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ABSTRACT 

We present high angular resohition (0'.'5 — 2'.'0) observations of the mm continuum and 
the i2cO(J=3-2), i3CO(J=3-2), i3CO(J=2-l), Ci80(J=2-l), HCN(J=3-2), HCO+(J=4-3) and 
HCO+(J=3-2) line emission in the circumnuclear disk (r<100 pc) of the proto-typical Seyfert type- 
2 galaxy NGC 1068, carried out with the Submillimeter Array. We further include in our analysis 
new ^^CO(J=1-0) and improved ^^C0(J=2-1) observations of NGC 1068 at high angular resolution 
(I'.'O — 2'.'0) and sensitivity, conducted with the IRAM Plateau de Bure Interferometer. Based on the 
complex dynamics of the molecular gas emission indicating non-circular motions in the central '^100 pc, 
we propose a scenario in which part of the molecular gas in the circumnuclear disk of NGC 1068 is 
radially blown outwards as a result of shocks. This shock scenario is further supported by quite warm 
(Tkin>200 K) and dense (n(H2)~10^ cm^^) gas constrained from the observed molecular line ratios. 
The HCN abundance in the circumnuclear disk is found to be [HCN]/[^^CO]wlO~^-^. This is slightly 
higher than the abundances derived for galactic and extragalactic starforming/starbursting regions. 
This results lends further support to X-ray enhanced HCN formation in the circumnuclear disk of 
NGC 1068, as suggested by earlier studies. The HCO+ abundance ([HCO+]/[i2Qo]salO-5) appears 
to be somewhat lower than that of galactic and extragalactic starforming/starbursting regions. When 
trying to fit the cm to mm continuum emission by different thermal and non-thermal processes, it 
appears that electron-scattered synchrotron emission yields the best results while thermal free-free 
emission seems to over-predict the mm continuum emission. 

Subject headings: Galaxies: individual: NGC 1068 - Galaxies: ISM - Galaxies: active - Galaxies: 
kinematics and dynamics - Galaxies: nuclei ~ Galaxies: Seyfert - Radio continuum: 
galaxies - Radio lines: galaxies - Submillimeter: galaxies 



1. INTRODUCTION 

Little is known about the effects of active processes 
in galaxies on the chemical and kinematic properties of 
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the surrounding molecular gas and vice versa, whether 
the activity is in form of an active galactic nucleus 
(AGN) or a starburst (SB) or both. Information on the 
characteristics of the molecular gas in the vicinity of the 
activity is essential to reveal the underlying physical 
processes because molecular gas constitutes a large 
fraction of the fuel for the central activity and thus helps 
to keep it alive over cosmologically relevant time scales. 
Also, the feedback of activity onto the surrounding 
molecular gas represents an important factor for the 
evolution of the activity (with respect to outflows or 
shocks for instance). The large diversity and oftentimes 
also simultaneity of the physical processes accompanying 
the different activity types certainly complicate any 
interpretation of the interaction between the activity 
and the molecular gas. These processes include large 
scale shocks, gas out- and inflow, other dynamical per- 
turbations, and strong radiation fields, such as through 
UV- or X-ray radiation, cosmic rays, or supcrnovae 
explosi ons (e.g.. iMartm et al.il201ll: iGarcia-Burillo et all 



Perez-Bca upuits et al.l 120091: iKrips et al 



20 1 ; 'S akamoto et a.l.| 20101: Papadopouloi 



2008; Matsushita et al."2007': 'Garcfa-Burillo et al .l 12007 
Martin et al. 2006; Usero ct_al, 2006; Saka moto et al.l 



2006t " (Matsus hita et a l.l l2005t , iFuente et al.l 120051: 
Meier fc Turner,. 2005; Usero et al.ll2004l ). ThuT a thor- 



ough study of the kinematics, excitation conditions and 
chemistry of the molecular gas close to AGN and SBs is 
essential for understanding the nature and evolution of 
these active environments. 
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High angular resolution observations of the ^^CO emis- 
sion, a reliable tracer of the global molecular gas reser- 
voir, are an important step to study the dynamics in ac- 
tive galaxies. ^^CO alone, however, can certainly not de- 
scribe the complexity of the molecular gas close to the ac- 
tivity processes (especially with respect to the chemistry 
and excitation conditions of the molecular gas). More- 
over, ^^CO has been found to be an unreliable tracer of 
dense molecular gas environments (at least its lower ro- 
tational transitions), in which most AGN or star forma- 
tion activity is supposed to take place (e.g., iKrips et"al] 
l2007HGao fc Solomon|[200l . Observations at high angu- 
lar resolution of different molecular tracers is thus a next 
logical step. Given recent upgrades of current (sub)mm 
interferometers, the detection and spatial resolution of 
weaker molecular lines becomes a feasible task. 

In this paper we present high angular resolution ob- 
servations of the (sub-)mm-continuum and ^^CO, ^'^CQ, 
C^^O, HGN and HCO+ line emission in the nearby 
Seyfert type-2 galaxy NGC 1068, conducted with the 
Submillimeter Array (SMA) and the IRAM Plateau de 
Burc Interferometer (PdBI). 

2. NGC 1068 

The nearby Seyfert type-2 galaxy NGC 1068 (see 
Table [1] for a general overview) has become not only 
the figurehead for the viewi ng angle unification the - 
ory for Seyfert galaxies (e.g. iKrolik fc Kallmaiil I1987D . 
as its centre most impressively exhibits the character- 
istics for harb oring an obscured type-1 Seyfert AGN 
(|Antonucci fc M iller 1985). It also advanced to a figure- 
head for the significantly different effects that AGN can 
have on the excitation conditions and chemistry of the 
surrounding molecu lar gas when compared to SB or qui- 
escent gala xies (e.g. Krips et al. 2008; Us ero et al.ll2004 
iKohno et a l. 2001; Sternberg ct al. 1994). The prototyp- 
ical nature of NGC 1068 is certainly in part due to its 
relatively small distance (12.6 Mpc) from us and strong 
continuum as well as line emission from X-ray to radio 
frequencies, making it hence an ideal target to study the 
accretion and feedback processes of its active nucleus at 
unprecedented detail. As a (fortunate) consequence, a 
wealth of information is already available for this source, 
the most revelant of which will be summarised in this 
section. 

The molecular gas in NGC 1068 is distributed in 
a starburst ring/spiral of ^-^3 kpc (~50") in diame- 
ter, a stellar bar of ~2 kpc (~30") in length and a 
cirumnuclear di sk/ring (CND) of ~2 00 pc (~3") in 
diameter (e.g., ISchinnerer et al.l l2000l and references 
therein). At the very center of the CND, pronounced 
II2O maser emissi on suggests a thin disk of a few parsec 
in dia meter fe.g.. iGallimore et al.ll200"ll iGreemhill et al.l 
|1996[) . A pronounced jet and counter-jet can be ob- 
served fro m cm to rnm w avelengths (e.g., iKrips et al.l 
[2006; Gall imore et al.|[2004l) that extends from the maser 
disk out to several kpc from the center. MIR observa- 
tions reveal hot and ionised gas that biconically follow 
the path of the radio jet ('e.g..lMiiller-Sanchez et al.ll2009l; 
Poncelet et al.''2008': 'T omono et al.l 120061; iGalliano et al.l 
2004; Bock ct al. 2000) a nd indicate the e xistence of a 
pc-scale warm dust torus (jJaffe et al.ll2004l ). Early high- 
angular resolution radio-continuum observations indi- 
cated an interaction of the radio-jet with the neighbour- 



ing insterstellar medium (ISM) due to the apparently 
disrupted structure of the (northern) jet (radio compo- 
nents NE & C in, e.g., IGallimore etal] ll996': 'Rov et al 
I1998D . However, in a later publication, IGallimore et al 
(2004) argue that the observed disturbed jet-structure 
(in NE) can as well be explained by variable outfiow 
speeds due to variable accretion (e.g., IGallimore et al.l 
I2OOII: ISiemiginowska fc Elvislll997 |) though an earlier col- 
lision at position C was not discarded 

The complex dynamics of the molecular gas as 
traced by the ^^C0(2-l) line emission within the 
CND were interpreted as a consequence of a warped 
disk (jSchinnerer et al.) l2000l ). However, more recent 
observations of the MIR rovibrational H2 emission 
start t o raise doubt s ab out this interpretation (e.g., 
I Miiller- Sanchez et al.l 12009) and alternatively suggest 
that the complex gas kinematics are due to a funneling 
of the gas toward the AGN along the jet (inner 60 pc) 
plus an expanding ring (on scales of r=100-150 pc), the 
latter having been already proposed a few years before 
by Gallia no fc A Uoin (2002). 

Additional fascinating characteristics of the CND of 
NGC 1068, besides the complex kinematic behaviour 
of its molecular gas, are its chemistry and excita- 
tion conditions which appear to significantly differ from 
starburst/star- forming environments. Traced by the 
"abnormal" line ratios of different molecules and tran- 
sitions, mostly by HGN, HCO+ and ^^CO, it has been 
suggested that the CND in NGC 10 68 harbors a giant 
X-ray- dominatcd-rcgion (XDR; e g., iKohno ct al.' loF" 
Usero et al. 2004; Tacc oni et al.l Il994; Sternberg ct al 
1994) . XDRs are defined in a similar way to the Photon- 



Dominated-R egions (PDRs) in s tarburst galaxies (such 
as M82; e.g., iFuente et all 120051 ) but are driven by X- 
ray rather than UV-radiation. High HCN-to-CO(J=l-0) 
(>1) and HCN-to-HCO+(J=l-0) (>1) ratios are found 
in the CND of NGC 1068, indicating enhanced HGN 
abundances there. The X-ray radiation of the AGN is 
thereby supposed to be the main driver for the enhance- 
ment of HCN. It can penetrate much deeper into the sur- 
rounding molecular gas than the UV-radiation in PDRs 
leading to a stimulated "hyper" -production of HCN. 
A multi-transition, multi-molecular line study of HCN 
and HCO+ condu cted with the IRAM 30m-telescope 
(jKrips et al.ll2008l) supports an increased abundance of 
HCN and/or increased kinetic temperatures. Both can 
equally explain the elevated HCN-to-^^CO(l-O) line ra- 
tios either by the aforementioned enhancement of the 
HCN abundance and/or a hypo-excitation of the low- J 
^^CO transitions (see also M51 and NGC 69 51 as exam- 
ples of incr e ased HCN/^^CO(J=1-0 ) ratio s : iKrips et al" 



19961) 



20091 120071; iMatsushita et all l2007l 12004 IKohno et al 



Recent SiO interferometri c observations of the CND 
in NGC 1068 carried out bv iGarci'a-Burillo et al.l (|2010f ) 
testify further to the complexity of the gas chemistry in 
this galaxy. The bright SiO emission in its CND suggests 
an enhanced abundance of this molecule which is inter- 
preted by the authors as closely related to (high- velocity) 
shocks. The shocks are believed to be a consequence of 
a jet-gas interaction. 

3. OBSERVATIONS 



Millimeter observations of NGC 1068 
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A summary of all (sub-)millimeter interferometric ob- 
servations is given in Table El where observing parame- 
ters as well as achieved rms noise and angular resolutions 
are listed. 

3.1. SMA 

For all SMA observations, unless otherwise 
stated, the phase reference centre has been set to 
aj20oo =02h42m40.70s and 5j2ooo =-00°00'47'.'9 which 
corre sponds to the radio po s ition of the active nucleus 
(e.g. iGallimore et~aLl 120041: iKrips et al.l I2006D . The 
SMA receivers have been tuned to the respective line 
using doppler-tracking on the systemic velocity of 
NGC 1068 of VLSR=1137 km s^^. Both the upper 
sideband (USB) and lower sideband (LSB) were used 
for the observations, yielding a bandwidth of 2 GHz 
each separated by 10 GHz. A spectral resolution of 
0.81 MHz was used for all observations, corresponding 
to 0.8 km s~^ at 1 mm. The 225 GHz zenith opacity T225 
was measured regularly throughout all observations at 
the nearby Caltech Submillimeter Observatory (CSO). 
The accuracy of the flux calibration for all tracks is 
estimated to be at a conservative level of ~20%. 

The SMA data have bee n reduced with SM A specific 
tasks in the MIR package ()Scoville et al.|[T993[) . Further 
image a nalysis has been conducte d with the GILDAS 
package (jGuilloteau fc Lucail2000[ ). 

3.1.1. HCN(J^3-2) emission 

We observed the HCN(J=3-2) line emission in 
NGC 1068 using the extended and very extended con- 
figurations with up to eight 6 m dishes in January, Oc- 
tober and November 2006. The 345 GHz receivers were 
tuned to the HCN(J=3-2) line (265.886 GHz at rest) 
in the LSB; the USB was used for continuum measure- 
ments. The weather conditions were good with opacities 
of T225=0. 05-0.1 in the January track (extended config- 
uration) and T225=0.13-0.22 in the October/November 
tracks (very extended configuration). We used 3C273, 
3C111 and/or 3C454.3 as bandpass and Uranus, Titan, 
and/or Neptune as flux calibrators^^. We observed two 
quasars (0235+164, 0238-084, 0339-017, and/or 0423- 
013) every ^15 minutes to calibrate the gains (amplitude 
and phase versus time). The data from all three tracks 
have been combined into one single data file, resulting in 
an rms noise of 24 mJy in 17 km s~^ wide velocity chan- 
nels. The synthesized beam is determined to be l'.'0x0'.'8 
at a position angle of PA=30° (natural weighting) and 
0'.'53x0'.'46 at PA=30° (robust weighting). 

3.1.2. HCO+ (.7^3-2) emission 

We carried out observations of the HCO''"(J=3-2) 
emission in NGC 1068 using seven antennas in extended 
configuration during November 2006. The 345 GHz re- 
ceivers were tuned to the HCO+(3-2) line (267.558 GHz 
at rest) in the LSB; the USB was used for continuum 
measurements. The weather conditions were good with 
opacities of r225=0.06-0.15. Bandpass calibration was 
performed on 3C273, Titan and Uranus, while absolute 
fluxes were determined using Titan. The gains have been 

We used the line-free USB respectively for Titan and Neptune 
to determine the absolute flux level as they are known to have broad 
HCN lines which could contaminate a flux calibration in the LSB. 



calibrated on 0423-013 and 0339-017. For this data set, 
we reach an rms noise of 33 mJy in 17 km s^^ wide ve- 
locity channels. The synthesized beam is determined to 
be I'.'OxO'.'S at PA=30° (natural weighting). 

3.1.3. ^'^CO(J^3-2) and HCO+ (J=4~3) emission 

The i2cO(J=3-2) emission of NGC 1068 was ob- 
served in extended conflguration using all eight anten- 
nas during September 2007. HCO+(J=4-3) was addi- 
tionally observed in a separate track in August 2007. 
The 345 GHz receivers were tuned to the ^^C0(J=3- 
2) line (345.796 GHz at rest) in the LSB such that the 
HCO+(J=4-3) (356.734 GHz at rest) stiU falls within 
the USB. The opacities ranged between r225=0.06-0.13. 
Bandpass calibration was performed on 3C454.3 and 
Uranus. Uranus was also used for flux calibration. Gains 
were determined using 0238+166 and 0423-013. An rms 
noise of 80 mJy is reached in 7 km s~^ wide velocity 
channels. The synthesized beam is determined to be 
I'.'OxO'.'S at PA=30° for natural weighting when also us- 
ing a uv-taper to better match the angular resolutions of 
the other SMA observations. The original (un-tapered) 
angular resolution amounts to 0.6x0.5 at PA=:30°. 

3.1.4. ^^C0(J^2'1) and C^^0(J^2-1) emission 

The i3C0(J=2-l) emission in NGC 1068 has been ob- 
served in extended conflguration using all eight anten- 
nas during January and February 2008. The 230 GHz 
receivers have been tuned to the ^■^C0(J=2-1) line 
(220.399 GHz at rest) in the USB such that the 
Ci^0(J=2-l) (219.560 GHz at rest), the HC3N(J=23- 
22) (209.230 GHz at rest) and H32a (210.502 GHz 
at rest) still fall within the LSB. However, only the 
13C0(J=2-1) and C^^0{i=2-l) line emission was de- 
tected. The opacities ranged between T225=0.1-0.2. 
Bandpass calibration was performed on 0423-013, 3C111 
and Titan, while gains were determined using 0339-017 
and 0423-013. Titan was further used as a flux calibra- 
tor. An rms noise of 12 mJy was reached in 17 km s~^ 
wide velocity channels. The synthesized beam is deter- 
mined to be l'.'0x0'.'9 at PA=30° for natural weighting 
when also using a uv-taper to better match the angular 
resolutions of the other SMA observations. 

3.1.5. ^'"^00(3=3-2) emission 

The i3C0(J=3-2) emission in NGC 1068 was observed 
in compact configuration using seven antennas in Oc- 
tober 2005. These observati ons were part of t he ob - 
serving campain presented by iHumphrevs et al.| ()2005D , 
which aimed to detect extragalactic H2O maser emis- 
sion at (sub-)millimeter wavelengths. The 345 GHz re- 
ceivers have been tuned to the H2O(10(2,9)-9(3,6)) maser 
fine (321.226 GHz at rest) in the LSB such that the 
i3cO(J=3-2) line (356.734 GHz at rest) was stifl lo- 
cated within the USB. The 225 GHz zenith opacity has 
remained stable around 0.05-0.06. Bandpass calibra- 
tion has been performed on 3C454.3, 3C111 and Uranus. 
Uranus has been also used as a flux calibrator. The gains 
were determined using 0234+285 and verified against 
0215+015 and 0420-014. An rms noise of 61 mJy is 
reached in 7 km s^^ wide velocity channels. The syn- 
thesized beam is determined to be 2'.'4x2'.'l at PA=28° 
(natural weighting). 
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3.2. IRAM PdBI 

3.2.1. ^'^C0(2-l) and ^^CO(J=1~0) emission 

Observations of the i2C0(2-l) emission in NGC 1068 
were carried out with the IRAM PdBI in February 
2003 using ah six antennas in A configuration. Si- 
muhaneously, we observed the ^^CO(J=1-0) using the 
3 mm PdBI receivers. The bandpass was cahbrated on 
NRAO150 and 0420-014 while phase and amplitude cal- 
ibration were performed on 0235-1-164 and 0238-084. A 
total bandwidth of 580 MHz with a spectral resolution 
of 1.25 MHz was used. We reach an RMS of ^7 mJy 
in 7 km s^^ wide channels (natural weighting) at 1 mm 
and of ^1.9 mJy in 14 km s^^ wide channels (natural 
weighting) at 3 mm. Applying natural weighting in the 
mapping process, beam sizes are derived to be l'.'0x0'.'6 
at PA=36° at 1 mm and 2'.'5xl'.'7 at PA=28° at 3 mm. 
However, to better match the SMA observations, we 
mapped the CO (2-1) data with a uv-taper giving an 
effective angular resolution of 1.0"x0.8" at PA=30°. As 
the uv-coverages between the SMA and PdBI data (for 
the high angular resolution data) are very similar, the 
usage of a simple uv-taper already provides the neces- 
sary accuracy to match the restoring beam of the PdBI 
observations with that of the SMA. 

4. RESULTS 

4.1. Continuum emission: from 850 fim to 1.4 mm 

The sub(mm) continuum emission at the wavelengths 
presented in this paper was derived from the line-free 
channels of the respective line observations (Table [2]), 
averaging emission from USB and LSB where possible 
(for the SMA data; the PdBI data were obtained from 
single-sideband observations). The continuum emission 
has further been merged between data sets with very sim- 
ilar observed frequencies (i.e., within ~10 GHz). Before 
averaging data from different sidebands and/or observa- 
tions, it has been carefully verified that the absolute flux, 
position and structure of the emission in the individual 
data sets are consistent with each other within the cali- 
brational uncertainties of 20% (flux) and 0.1" (position) 
in order to reduce systematic errors. 

The continuum emission is clearly detected at all wave- 
lengths at a >5cr level. To obtain accurate fluxes, posi- 
tions and sizes, elliptical Gaussians were fltted to the 
data in the uv-plane, except for the uniformly weighted 
map of the 1.0 mm continuum emission for which a cir- 
cular Gaussian was fitted given its apparently unresolved 
nature. The results of these fits are listed in Table [H 

The continuum emission at 1.0 mm (NA), 1.3 mm and 
1.4 mm appear to be consistent with each other in terms 
of their flux, position and struc ture (see Table [3] and 
Fig. m see also lKrips et al.ll2006D . All show peak fluxes 
of around 15-19 mJy/beam and spatially integrated flux 
densities of 22-28 mJy, indicating extended emission. 
Their positions, although self-consistent, are slightly to 
the North (^0.2" ) of the radio position o f the AGN (com- 
ponent SI from iGallimore et al.l l2004l marked with a 
white cross in Fig. [1]) and that of the uniformly weighted 
1.0 mm continuum emission (Fig. [T}d; white contours). 
The shift between the mm and cm data is larger than 
the positional uncertainty of 0.1" and thus assumed to 
be real. NGC 1068 is known to have a pronounced radio 
(and mm-) jet in a North-East-to-South-West direction. 



of which the North-Eastern part exhibits the st ronger 
emission (e.g., IGallimore et~al] 12004 iKrips et al.l [20061. 
Despite the steep spectral index of the synchrotron emis- 
sion of the jet, the extended (i.e., >1") emission from 
both the jet and counter-jet are still visib le at 3 mm (e.g., 
IKrips et al.ll2006l : ISchinnerer et al.ll2000f ) but are signifi- 
cantly fainter or undetected at shorter wavelengths (i.e., 
<1 mm). 

The continuum emission at 1.0 mm (NA), 1.3 mm and 
1.4 mm is a blend of emission associated with the (North- 
East) radio jet and the AGN itself (SI in lGallimore et all 
l2004f ) due to the 'lower' angular resolution of ^^1". Due 
to the higher angular resolution, the 1.0 mm contin- 
uum emission of the jet in the uniformly weighted map 
(Fig. [TId) is almost entirely resolved out leaving behind 
only the more compact emission of the AGN. Thus, the 
centroid of the emission at 1.0 mm (NA), 1.3 mm and 
1.4 mm will naturally be shifted towards the North, while 
the 1.0 mm (UN) continuum emission should reveal the 
actual position of the AGN (or at least the base of the 
jet). 

Going to even shorter wavelengths of 850 /im, it ap- 
pears that not only the continuum flux increases again, 
but also its position seems to be now consistent with 
the AGN, independent of the weighting (i.e., synthesized 
beam) used for mapping/cleaning and unlike the 1.0 mm 
(NA), 1.3mm, and 1.4 mm continuum emission. The 
latter may indicate that the emission from t he radio jet 
is ne gligeable at 850 /zm (see also Fig. 3 in IKrips et al.l 
l2006l ) and the AGN (i.e., the SI component) dominates. 
The increased flux at 850 /im, which appears to be larger 
by almost a factor of 2 compared to the 1.0 mm-1.4 mm 
emission, shows that thermal dust emission already plays 
a signiflcant role at 850 /im (see Section lO]) . Also, the 
size and shape of the continuum emission appear to have 
changed compared to that at longer wavelengths. The 
PA of the 850/im emission of ^^90° is signiflcantly differ- 
ent from that ('^30°) of the 1.0 mm-1.4 mm emission. 
Moreover, the 850 /xm emission appears to be extended 
(Fig. [TJ;), in contrast to the 1.0 mm-1.4 mm emission, 
which seems to be extended only in the jet component 
but not in the 'left-over' AGN component in the uni- 
formly weighted 1.0 mm map (Fig. [T]d). The uniformly 
weighted 850 /im continuum emission appears to be also 
resolved (white contours in Fig. [TJ;, compare also peak 
flux with total flux density in Table [3]) . 

4.2. Line emission 

4.2.1. General Characteristics and Distribution of ttie 
Molecular Gas 

The continuum emission has been subtracted from all 
line data in the uv plane to avoid any contamination of 
the line by continuum emission even if in some cases the 
continuum emission does not exceed the noise level in 
the individual channel maps (see Table [2] & [3]) . In or- 
der to reduce systematic effects due to spatial filtering, 
we used a slight uv-taper (giving some more weight to 
the shorter baselines) to map and clean all line emission 
data with the same synthesized beam, except for the high 
angular resolution ('^0.4-0.5") of the uniformly weighted 
HCN(J=3-2) map (shown additionally in Fig. [2^ using 
white contours) and the low angular resolution (^2") 
of the i2co(J=l-0), 13CO(J=1-0) and ^^€0(3=3-2) 
maps. 
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Fig. [2] shows the velocity integrated intensity maps 
of the molecular line emission from HCN(J=3- 
2), HCO+(J=3-2), HCO+(J=4-3), i2co(J=2-l), 
13C0(J=2-1), "CO(J=3-2), Ci80(J=2-l), i2cO(J=l- 

0) and "CO(J=1-0); i2cO(J=3-2) is plotted in grey 
scale in all images to facilitate a comparison. All these 
molecules have been clearly detected above the 5a level 
(except ^'^CO(J=1-0)). The emission in all lines reveals 
a pronounced peak on the stronger eastern knot and in 
the stronger lines also a weaker peak on the western knot, 
both of which are already known from previous ^^CO 
observations (e.g.. ISchinnerer et al.ll2000l ). An elliptical 
Gaussian has been fitted to the uv-data for all lines in 
order to obtain the position, peak- and spatially inte- 
grated flux of the emission in the two knots. The re- 
sults of the fits are given in Table ID The position of the 
emission in the eastern and western knot is very simi- 
lar in all observed lines, excluding the ^'^C0(J=2-1) and 
C^^0(J=2-1) emission which seem to peak closer to the 
AGN. 

The spectrum of the spatially integrated emission (over 
the central '--^4") of each line is plotted in Fig. [3l We 
also show the ^ ^ CO(J =1-0) line emission taken from 
ISchinnerer et al.l ()2000D for consistency. While the veloc- 
ity integrated line emission seems to be very similar in its 
shape and position for most lines, the line profiles vary 
significantly from each other. While for the dense gas 
tracers (HCN, HCO"'") a single Gaussian fit is sufficient 
to reproduce the line profiles, the CO lines need a dual, 
triple or quadrupole Gaussian fit. However, to simplify a 
comparison, the results given in Table[S]represent a single 
Gaussian fit to all lines. The line centers are roughly con- 
sistent with each other, differing by less than 20 kms s~^. 
Excluding the ^^CO(3=2-l) and 0^^0(3=2-1) emission, 
the line widths also agree with each other within the un- 
certainties. Except for the i2C0(J=2-l) and "C0(J=2- 

1) line emission for which roughly half of the emission 
seems to be resolved out, the interferometric observations 
have captured most of the emission measured with sin- 
gle dish obser vations (T able [St . Please note that for the 
12CO(J=1-0) (ISchinnerer et al. 2000) and "CO(J=1-0) 
emission, the single dish fluxes are much higher than 
the interferometric ones because they contain signiflcant 
emission from the star-forming ring/spiral-arms and the 
bar. 

4.2.2. Dynamical Characteristics of the Molecular Gas 

The kinematic behaviour of the different molecules is 
presented in detail in Fig. 2] to [TT] To better under- 
stand the puzzling complexity of the different proflles of 
the various molecular lines and test whether it is due 
to dynamical effects, we spatially split the spectra by 
deriving the spectrum of the western and eastern knot 
separately (Fig. g]). Please note that the i3CO(J=l-0) 
and C^^0(J=2-1) line emission were discarded because 
of their insufficent sensitivity and/or lack of emission in 
the western knot while 12CO(J=1-0) and ^^CO{3=3-2) 
are not included because of their insufficient angular reso- 
lution. The iso- velocity maps (Fig.E]) of the ^^CO, ^^CO, 
HCN and HCO"'" line emission clearly show a dynamical 
structure that seems to be dominated by standard disk 
rotation with a blueshifted eastern knot and a redshifted 
western knot. If disk rotation were the only underlying 
kinematics, one would expect to flnd a simple blueshifted 



peak at the eastern knot and a redshifted peak at the 
western knot. Although disk rotation is observed. Fig. |4] 
shows kinematic features signiflcantly differing from sim- 
ple rotation. Instead, the blueshifted eastern knot also 
exhibits redshifted emission and the redshifted western 
knot blueshifted emission. These 'wings' appear to be 
present in all three ^^CO emission at a high signiflcance 
level as well as in the ^^CO, HCN and HCO+ emission 
but, given the lower signal-to-noise ratio (SNR) for these 
lines, not as pronounced as for ^^CO. At this point, it 
should be emphasized that in such a case the moment one 
map can be very misleading as it derives only the domi- 
nant kinematic structure and might overlook more com- 
plex underlying kinematics. Integrating (in velocity) the 
red- and blue-shifted parts of the line spectrum (Figl6]) as 
well as analysing the channel maps (Fig. [7]&[51) might be 
the more appropriate approach. Fig. [5] indicates a more 
complex distribution than expected from simple disk ro- 
tation. We flnd blueshifted emission spatially coincident 
with redshifted emission and vice versa; this seems to be 
most pronounced for the ^^CO, ^^CO and HCN emis- 
sion. By looking at the channel maps of the ^^C0(J=2- 
1) and ^^C0(J=3-2) emission (which have the highest 
SNR), the red- on blueshifted and blue- on redshifted 
emission is not only at low-velocities but also at higher 
velocities (which is especially visible in the ^2CO(J=2~ 
1) emission; see channels <— 70km/s and >-|-100 km/s in 
Fig. [7]). The same is true for the ^^CO(J=3-2) emission 
(Fig. [51 although less pronounced, especially for veloc- 
ities >80 km s""'^ for which no emission can be found 
anymore as opposed to ^^C0(J=2-1)). We flnd a be- 
haviour of the HCN(J=3~2) emission similar to that of 
the i2C0(J=2-l) and ^^€0(3=3-2) emission though on 
a much lower significance level. The ^^C0{3=2-1) emis- 
sion seems to indicate, however, a different behaviour 
(see Fig. [S]). Instead of being distributed in a 'ring'- 
like manner, the emission appears to be elongated more 
in a South- West to North-East direction (see especially 
channels maps between -1-50 km s^^ to —20 km s^^). 
However, given the low sensitivity level, this structure 
has to be treated with caution and needs confirmation 
by either higher sensitivity observations or other molec- 
ular lines such as SiO. Indeed, the SiO emission seems 
to indicate a similar behaviour as discussed in separate 
papers (see iGarcfa-Burillo et al.l[2008l[2010l) . 

The position-velocity diagrams of the ^^CO and HCN 
emission, taken at different Position Angles (PA) in steps 
of 30° across the CND, are shown in Fig. [TUl The grey 
scale denotes the ^^C0(J=2-1) emission for better com- 
parison. Overall, the kinematic structures in the different 
lines strongly resemble each other. Also, the overlap of 
the red- on blueshifted emission can be seen quite well 
in the position- velocity diagrams (see especially PA=60- 
120° in Fig. [TU)) . strongly indicating pronounced non- 
circular motions in the CND of NGC 1068. 

In order to quantify and parametrise the observed 
comp lex kinematics, we f ollow the approach used by 
iHeckman et al.l (|1989[ ) and iBaum et al'l (|1992i ). We de- 
termine three kinematic parameters from the slits taken 
at the different position angles used in Fig. [10] for the 
^^CO emission. These parameters are: 1.) the av- 
erage linc-of-sight velocity dispersion cr, determined as 
0.426xFWHM along each slit, 2.) the "rotational" veloc- 
ities A, determined from the difference between the aver- 
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age velocities on either side of the nucleus along each slit, 
and 3.) the rms variation of the velocity e for each point 

along the slit, defined as e = \J jfT,fLi{vi — Uav)^: where 

N is the number of points along the slit, Vi the intensity 
weighted velocity for point i and the Wav the intensity 
weighted average velocity along the slit. A comparison 
of the different parameters with each others, especially 
the ratios ^ and ^ allows to classify the dynamics into 
three differen t groups (for a more detailed explanation 
we refer to Baum et alJ (Il992[ )): 

ROTATORS: A/cr > 1, A/e > 1 

CALM NON-ROTATORS: A/a < 1, A/e - 1 

VIOLENT NON-ROTATORS: A/a < 1, A/e < 1 

The results of this kinematical parametrisation for the 
^^C0(J=2-1) are shown in Fig. [11] highlighting our pre- 
vious findings (we find very similar values when taking 
the ^^CO(J=3-2) emission). For most position angles, 
we find that the ratios are most consistent with calm non- 
rotation with one exception (PA=60°) being located in 
the area of the violent non-rotators. This strongly em- 
phasizes the fact that although there is an underlying 
dominating disk rotation, the dynamics of the molecu- 
lar gas in the CND of NGC 1068 is significantly dis- 
turbed by a non-rotational process, most significantly for 
the North, North-Eastern part of the CND (i.e., along 
PA=0-90°). We also attempted to determine the rota- 
tional velocities directly from azimuthally averaging the 
velocities within the CND and then subsequently fitting 
rotational curves to the data using rotcurv in gipsy. We 
thereby assumed different scenarios ranging from pure 
disk rotation to adding a radial dependency. However, 
no physically meaningful result could be obtained. We 
mostly find positive velocities that decrease with radius 
which is inconsistent with simpl e disk rotation (see also 
the dynamical analysis done in ISchinnerer et al.l l2000l) 
and necessitates the inclusion of a bulge or central mass 
component whether in form of nuclear star cluster or a 
massive black hole in addition to the disk. This certainly 
further emphasizes the complexity of the molecular gas 
dynamics in the CND of NGC 1068. 

4.2.3. Molecular Line Ratios 

In order to constrain the excitation conditions and 
chemistry of the molecular gas, we derive the line ra- 
tios for the different molecules and transitions in sev- 
eral ways, by accounting for the different angular res- 
olutions (especially with respect to the ^^CO(J=1-0), 
13CO(J=1-0) and i3C0(J=3-2) emission). Before de- 
termining any line ratio, the line emission was brought to 
the same (lower) angular resolution by using a uv-taper; 
this seems appropriate for most of the lines as we recov- 
ered most of the emission with the interferometric obser- 
vations. Also, we compute line ratios only for emission 
coming from the samwe spatial regions (see Figs lT2lfT4|) . 
Fig. [12] shows the velocity integrated line ratios between 
various combinations of the molecular lines. Separating 
spatially the eastern and western knot we derive spatially 
averaged line ratios from Fig. [T^] which are listed in Ta- 
ble [5] Please note that the values in Table [Bj might vary 
from those derived from Tables [4] and [5j However, the 



differences can be easily explained by the different se- 
quence of averaging (i.e., first in space, then in velocity 
versus first in velocity, then in space) for Table [SJ the us- 
age of an elliptical Gaussian fit^^ for Table [4] as opposed 
to spatially averaging without a fit as done for Table [6l 
and the different spatial resolutions of the line emission 
in Tables [1& [5] 

We identify some spatial variance of the different line 
ratios (mostly a factor ~2-3) which seem to be most pro- 
nounced in the ^^CO and ^'^CO hue ratios (Fig.fT^.d.o). 
The HCN and HCO+ line ratios seem to be more con- 
stant over the CND than ^^CO and "CO. The higher 
values are found closer to the position of the AGN for 
most maps. 

In order to investigate whether there might be a veloc- 
ity (and spatial) dependence on the line ratios, we deter- 
mined the line ratio channel maps for the two strongest 
^^CO transitions (J=2-l and J=3-2) as function of ve- 
locity in Fig. [T31 We find somewhat higher (i.e., fac- 
tor of >2) i2co(J=3-2)-to-i2cO(J=2-l) ratios close 
to the AGN at velocities around the systemic velocity 
but also on the eastern knot at high negative velocities 
(< — 130 km s~^). Both knots seem to show (more or 
less) the same velocity and spatial behaviour as can be 
seen in Fig. [141 This plot shows the spatially averaged 
i2CO(J==3-2)-to-i2CO(J=2-l) line ratios for the east- 
ern and western knot as function of velocity. The two 
curves follow each other nicely except for velocities be- 
tween — 120 to —140 km s^^ for which the eastern knot 
exposes higher values (by a factor of 2). The error bars 
denote the variance of each averaged value which in most 
cases indicates a variation by a factor of 1.5. 

5. DISCUSSION 

5.1. Spectral Energy Distribution of the Continuum 
Emission 

The nature of the continuum emission (from IR over 
sub-mm to cm wavelengths) represents a highly debated 
and complicated matter for NGC 1068, recently gain- 
ing a revival by newly published VLTI/MIDI (i.e., IR) 
and radio data (e.g., iHonig et al.l l2008t ICotton et al.l 
l2008f ). As mentioned in the previous section, the ra- 
dio continuum emission splits up into several compo- 
nents, a jet plus counter-jet and a core component (SI) 
associated with the AGN itself. While the emission 
from the jet is certainly pure non-thermal synchrotron 
emission, a s supported by its st e ep continuum spec- 
trum (e.g., iGallimore et~aLl 120041: ICotton et all [2008), 
the natu re of the emissi o n fro m SI is highly contro- 
versial. IGallimore et al.l (|2004D already rule out syn- 
chrotron emission as origin for the continuum spectrum 
of SI and discuss electron-scattered synchrotron emis- 
sion as well as thermal fre e -free absorption as alter- 
natives. While iKrips et al.l (|2006[ ) present arguments 
for electron-scattered synchrotron emissio n based on a 
turnover seen between cm and mm-data, ICotton et al.l 
(2008) rather support the thermal free-free absorption 
model. A highly complicating factor in this discussion is 
certainly the mismatch in angular resolution between the 

Although this a reasonable first order fit, the line emission is 
certainly not exactly of an elliptical Gaussian shape so that some of 
the emission is not well reproduced by fitting an elliptical Gaussian 
to the velocity integrated maps. 
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cm, mm and IR data. As discussed in lKrips et al.l (|2006[ ) 
and this paper, the mm-continuum emission is contami- 
nated by emission from the jet at angular resolutions of 
>1", introducing large uncertainties in the estimate of 
Si's flux (compar e iHonig et al.ll2008HCotton et alJl2008l : 
iKrips et"al1l2006( ) due to the lack of angular resolution. 
However, the previous estimate of the 1.3 mm continuum 
flux of (10±4) mJy in lKrips et al] (|2006[) . translating to 
(9±4) mJy at 1.0 mm, is very similar to the measured 
1.0 mm continuum flux of (13±2) mJy from our high- 
angular resolution SMA observations. Given the unre- 
solved nature of the latter, jet emission does not seem 
to be signific a.nt anymore a t this angular resolution (see 
also Fig. 1 in lCotton et al.| [2008) although the obtained 
angular resolution is still an order of magnitude larger 
than that at cm wavelengths. 

Taking also the new 850 ^m (UN) continuum flux mea- 
surements into account, we recomputed the spectral en- 
ergy distribution (S EP) and replotted the m odels used by 
IHonig et al.l (|2008l ) and lKrips eraD ()2006f l. We thereby 
base our graphs on the formula and p arameters specified 
in Equations (2)- (5) and Table 2&3 in'Honig et al. (2008) 
and Equations (l)-(2) and Fig. 3 in iKrips et al.. (,200& ). 
The results are shown in Fig. IT5k-c. We marked all data 
points with a circle for which the obtained angular reso- 
lution of the observations did not exceed 1"^^. We also 
fitted a two-temperature grey body to the IR data, in or- 
der to estimate the contribution of thermal dust emission 
to the sub-mm continuum emission. Although this grey- 
body fit is certainly not as s ophisticated as the clumpy 
torus model used in IHonig e t al. ( 2008) , it represents a 
reasonable approximation as demonstrated by the good 
match to the IR data points. 

As can be seen in Fig. [T5b and c, the models used 
bv IHonig et al.l (1200 81 significantly overestimate the ob- 
served 1 mm (UN) flux by a factor of 2-3, although 
they reproduce correctly the 850 ^m one. It seems 
that the electron-scattered synchrotron emission model 
in Fig. 115b and the thermal free-free emission model in 
Fig. [T5b both need the extra contribution from the ther- 
mal dust emission to correctly reproduce the 850 /^m 
(UN) flux, while the synchrotron model in Fig. ITSb 
does not require it. Thus, it appears to be indeed very 
likely that the continuum emission is dominated by ther- 
mal dust emission starting at wavelengths <850/im, as 
posited above. 

Based on the 1 mm (UN) flux, it seems that the model 
best reproducing the SED at cm and mm wavelengths 
is the electron-scattered synchrotron emission; the ther- 
mal free-free absorption seems to overpredict the 1 mm 
flux. However, observations of the continuum emission in 
NGC 1068 have to be conducted at similarly high angu- 
lar resolutions (^0.5"), in order to dispel all remaining 
doubts, although the new mm observations presented in 
this paper are already a step in the right direction. 

5.2. Dynamics of the Molecular Gas 

In previous studies, the complex kinematic behaviour 
of the molecular gas has been thought to be a conse- 
quence of a warped disk. The warped disk has been 

^® Please note, that this is true for all data points except that at 
3 mm. The encircled SI data point has been estimated at 3 mm, 
not observed. We added this data point for consistency reasons 
only. 



modelled with a tilted ring model (e.g.. ISchinnerer et al.l 
I2OOOI ). However, the spatial overlap between the red- and 
blueshifted emission (i.e., the existence of highly non- 
circular motions) cannot be reproduced by these tilted 
ring models because they are based on circular motions 
and thus cannot account for non-circular motions of the 
gas (within the plane). 

Even though we cannot rule out a warped disk sce- 
nario in which part of the gas could be trapped in el- 
liptical orbits producing the non-circular motions, we 
want to propose an alternative approach, following re- 
cent findings on the H2(l-0) S(l) emission at scales of 
100-150pc bylMiiller-Sanchez et al.' ('200911) and the model 
proposed by Galliano & Alloin ( 2002) . The nature of the 
dynamics displayed in Fig. IH [Si |6l [7] and [H could also 
be explained by the following scenario: a rotating disk 
plus an outflow of the disk gas due to shocks and/or a 
CND-jet interaction. This hypothesis seems to further 
gain support when considering b esides the H2 1-0 S(l) 
map (jMiill er-Sanchcz et al.' '2009'. see their Fig. 4), also 
the 12/zm map (Bock_ct al. 2 00(1, see their Fig. 4) and 
the 5cm radio-continuum map (iGalhmore et al.ll2004 see 
their Fig. 1). The H2 1-0 S(l) and 12/zm emission follow 
nicely that of the radio jet in the inner 1" (North/North- 
East direction) which seems to interact with emission 
from the molecuar gas in the CND at ±1-2" in the North- 
ern part (see next Sub-Section); both the blue and red- 
shifted components associated with the non-circular mo- 
tions to the East and West of where the jet goes through 
or lies in front o f the CND (see also the case of M51; 
Mats ushita et al.l l2007. 2004). As the jet shows a biconi- 
cal structure with a change in direction close to the CND, 
it is not unreasona ble to believe that par t of it indeed hits 
the CND (see also lKraemer et aLlHOQi ). Such an inter- 
action could easily produce a shock through/along the 
CND (at least in the northern part, i.e. the bridge be- 
tween the eastern and western knot) feeding the assump- 
tion that some part of the molecular gas in the ring/disk 
might be blown outwards. 

Other causes for expanding/shocked gas include hy- 
pernovae explosions, stellar winds from a super stellar 
cluster as suspected in some n earby starburst galax- 
ies (such as NGC 253 or M82, ISakamoto eFall 120061 : 
iMatsushita et al.l l2000l ) or cloud-cloud collision within 
the CND (i.e., within the inner Lindblad resonance; see 
iGarcfa-Burillo et al.ll2010l ). However, they seem rather 
unlikely since the CND of NGC 1068 does not show any 
signs of starburst activity and the expansion seems to 
be too "ordered", i.e., too symmetric, to be caused by 
highly chaotic cloud-cloud collisions. Also, we cannot 
entirely rule out that the dynamics we see in the molec- 
ular gas is due to inflowing rather than outflowing gas, 
especially since indications of inflowing gas along the jet 
within the central ~70 pc (i.e., on scales smaller than the 
CN D) have been already prese nted in previous studies 
(see lMiiller-Sanchez et al.1l2009D . However, based on the 
appearance of the molecular gas within the disk (ring- 
like structure with an apparent void of gas in the inner 
part) lets us favor the jet-gas interaction rather than an 
inflow scenario (on scales of 100-150 pc). An outflow 
on scales of 100-150 pc is not necessarily in contradic- 
tio n with an inflow scen a rio on scales <70 pc proposed 
bv lMiiller-Sanchez et al] ()2009D . The jet- gas interaction 
could equally drag gas outwards on scales larger than 
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100 pc but trigger an inflow at smaller scales, depending 
on the type of interaction between the jet and the molec- 
ular clouds. However, higher-angular resolution obser- 
vations at sub-arcsecond angular resolution (as possible 
with ALMA) will certainly help to distinguish between 
the different scenarios. 

Based on this alternative approach, we tried to repro- 
duce the molecular gas dynamics with a very simplistic 
model that includes a dominant (Keplerian) rotating disk 
plus an outflow of some of the disk gas. We thereby as- 
sume a velocity gradient of Avo=200 km s~^, a radius of 
~1.2" and an inclination of ~60° for the rotating disk. 
We additionally add a slight ellipticity of 5% and assym- 
metry ratio between the eastern and western part of 1:0.7 
to the model. The outflow is approximated by a slightly 
expanding elliptical ring. We assume the same elliptic- 
ity of 5% as before and that the expansion starts at the 
inner radius of the disk. The expansion rate (H) is cho- 
sen to be 200 km s~^ per 1" or 67 pc equivalently (i.e., 
H=3 km/s/pc); this value is similar to w hat has been es - 
timated for M51 (2.2 km/s/pc; Matsus hita et al.ll2007l ). 
We also introduced a slight asymmetry ratio of 1:0.9 be- 
tween the blue- and redshifted emission. We further as- 
sume that up to ~30% of the disk gas is expanding. Most 
of these values that were chosen for the model parameters 
were almost instantly obvious from the observations, es- 
pecially the velocities and radii. We hence used them as 
starting values and scanned then through a reasonable 
parameter space for the optimal combination of input 
values that matches best the observed maps. However, 
given the larger number of parameters used in this model 
and hence the many degrees of freedom, we did not ac- 
tually conduct a true fit to the data but rather a "fit" 
by eye. The results of this so found "best-fit" model are 
displayed and compared to the ^^C0(J=2-1) emission 
in Fig. [TBll^ Indeed, the major dynamical character- 
sties of the molecular gas emission can be reproduced 
by this simplistic model supporting the hypothesis of an 
additional gas outflow. Also, using the same kinematic 
parametrisation for the model as used for the ^^CO emis- 
sion (see Fig. llip . we find very similar ratios between the 
observed and simulated emission. We have to stress that 
this suggestion does not exclude a warped CND but it is 
not needed for this model. 

It is interesting to note as further support for our ap- 
proach that the distribution of the HCN(J=3-2) emission 
matches almost exactly that of H2 1-0 S(l) in the CND, 
even better than the ^^C0( J=2-1) emission (see Fig. 1 in 
I Miiller- Sanchez et al.ll2009l) . Similar to the H2 map that 
indicates the brightest emission toward the North, also 
most of the HCN emission is found toward the northern 
part of the CND at which most of the potentially shocked 
gas would lie. Thus, one would expect the densest (and 
hottest) part of the gas in that area (see next section). 

5.3. Excitation conditions of the gas 

The line ratios derived from the interfcrometric maps, 
especially for the HCN, HCO+ and ^^CO emission 
(Fig. [T2] and Table 0, are consistent wi th previous 
findin gs from single-dish observations (e.g., iKrips et"al] 
|2008() . They support a picture in which the molecu- 
lar gas in the CND is relatively dense (n(H2)<10'*'^) 
and warm (Tk>40 K) with potentially higher than 
normal (i.e., in galactic giant molecular clouds) 



HCN abundances (Z(HCN)=[HCN]/[H2]) of Z(HC N)=1- 
50 xZgaiacMc(HCN)(Zgaia4tic (HCN) =2x10-8; e.g.. l!r^ 
119871 ). iKrips et all (|2008[) argue that the high HCN-to- 
12CO(J=1-0) and HCN-to-i3co(J=l-0) in NGC 1068 
can be explained by either higher t han normal HCN 
abundances due to an XDR (see also lUsero et al.l 120041 : 
ISternberg et al.l 119941) and/or higher gas temperatures 
leading to hypo-excited CO(J=1-0) emission. The lat- 
ter is supported by decreasing HCN-to-CO line ratios 
with increasing rotational number J. However, strong 
constraints on the kinetic temperatures could not be set 
based on the single-dish observations alone. Also, most 
single-dish observations are unable to unambigiously dis- 
tinguish between the molecular gas emission in the center 
and that in the star-forming spiral arms, complicating 
any interpretation of the data. Furthermore, most of the 
interfcrometric data previously published mostly focus 
on ^^CO at moderate angular resolution and only two of 
its transitions. 

The new interferometric maps, obtained for several 
transitions and molecules at sufficiently high angular res- 
olution, overcome some of the short-comings of previous 
observations / analyses. 

In the following we will discuss results from simulations 
of the excitation conditions of the molecular gas carried 
ou t with the radiative transfer code RADEX developed 
bv IVan der Tak et all (|2007[) . Please note that we did 
not find significant difference when using the LVG code 
in MIRIAD or the RADEX code. Given simplified simu- 
lations with RADEX, we decided to use RADEX in this 
paper as opposed to Krips et al. (2008) in which an LVG 
code was used. 

RADEX offers three different possibilities for the es- 
cape probability method: a) a uniform sphere, b) an ex- 
panding sphere (Large- Velocity-Gradient, LVG), and c) 
a plane parallel slab (shock) . We used all three methods 
but did not find significant differences for our data be- 
tween them. Thus, in order to keep the interpretation as 
simple as possible, we will discuss the results with respect 
to a uniform sphere in the following. 

We carried out simulations with RADEX for each 
molecule using the following grid of parameters (dimen- 
sion: 51x51x51): 

• Kinetic Temperature: Tkin= 1 - 500 K 

• Gas density: n(H2)=103 - lO"^ cm^^ 

• Column density: N(i2C0) = 10" - lO^i cm^^ 

We define the abundance ratio be- 
tween one molecule (MOLl) and another 
(M0L2) as X|;^gLi=Z(M0Ll)/Z(M0L2) with 
Z(M0L1) = [M0L1]/[H2]. For the different molecu- 
lar abundance ratios we assume the following ranges: 

• Abundance ratios: 

Xiagg ^ 6-1000 
XHCN^ ^ 1_500 

The ranges were chosen such that they span the 
abundance ratios found in different (galactic and ex- 
tragalactic) environments. In the Milky Way typical 
values of the ^^C/^'^C abundance ratio are ~20 which 
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increase to values of 80-100 in the outer parts of the 
galax y (e.g. IWilson fc RoodI 11994 IWilson fc Matteuccil 
|1992|) . Nearby starburst galaxies and ULIRGs show 
somewhat higher ^^C/^^C abundance ratios of >30 than 
found in the Galactic Center (e.g. NGC 253 M82 
IC 342, NGC 4945, NGC 6240 see |Greve et all IMil 
Henkel et al.lll998l[l99lll993al: iHenkel fc Mauersberger 



1993bD with Arp 220 being the exception. iGreve et al' 
(|2009[ ) find a very low abundance ratio of only 8 for 
this galaxy. Furthermore, some high redshift galax- 
ies seem to exhibit also rather high ^^C/^^G values of 
>30 as determined for the ISM in the graviational lens 
of P KS 1830-211 (jMuller et all l2006l ) or the Clover- 
leaf (jHenkel et al.l |2010D . Values of X^co+ 
equally large scatter ranging from the order of unity in 
starforming reg ions in the Milky Way (such as Orion 
and SgrB2; see lBlake et al.llT987D and nearby starburst 
galaxies/ULIRGs (such as M82, NGC 214 6, NGC 253, 
NGC 4945, NGC 6240 and Arp 220; see iNavIor et al 



2010t IGreve et al.l 120091: iKrips et al.l 120081: IWang et al 



20041 ) to >10 in nearby Seyfert galaxies (^e.g.. lKrips et al 
2008). 

For the simulations, we will concentrate on the region 
of the CND that contains all molecules. This region cor- 
responds to the bridge between the eastern and western 
knot, i.e., the northern part of the CND, and has a size 
of roughly -2-3" (~120-180 pc). This corresponds to 
roughly up to 50% of the molecular gas in the CND. If 
jet interaction indeed plays a role in NGC 1068, this will 
be the region most obviously affected. 



5.3.1. 



'CO & "CO emission 



We conducted a x^-test on the RADEX grid by using 
the observed line ratios for the ^^CO and ^^CO line emis- 
sion. Fig. |5T] shows the parameters for the best x^-test 
for four exemplary abundance ratios (10,26,52,110) from 
the aforementioned range; abundance ratios at the lower 
and higher end of the range show somewhat higher x^ 
values. 

The middle panel shows the lowest x^ found in each 
range of column densities and abundance ratios. The 
lower and upper panels show the respective lower and 
upper limits for the column density for which still a rea- 
sonably low x^ was found to indicate the spread in col- 
umn densities. 

The availability of the three lowest transitions for ^^CO 
and ^'^CO allows us to set tight constraints on Tkin, 
n(H2) and N(i2C0). The observed CO line ratios most 
impressively restrict the kinetic temperatures to values 
well above 200 K. This strengthens previous indications 
of warm/ho t ( Tkin > 50 K) mol ecular gas in the CND of 
NGC 1068 (iKamenetzkyetlll (Milt iKrips et al.l 120081: 
iMatsushita et al.l 119981: ISternberg et al.l ll994D~but be- 
ing rn uch higher t han the value f ound by lTacconi et al.l 
() 19941 ). However, iTacconi et al.l (|1994[ ) base their sim- 
ulations on single-dish data which cannot distinguish 
well between emission from the starforming ring/spiral 
arms and the CND. Especially the ^^CO(J=1-0) emis- 
sion might be overestimated for the CND which clearly 
is hardly detected in the interferometric map despite the 
high sensitivity of the observations. The lowest x^ is 
actually found for kinetic temperatures around 450 K, 
which seems to be fairly high. Although we used a 
one-gas component model due to the lack of sufficient 



observational constraints (i.e., higher-J CO transitions 
with Jupper>4), we do not expect all of the gas to be 
at these high kinetic temperatures, but rather a frac- 
tion of the gas in the northern part of the bridge. It 
seems likely that by using a two-temperature model for 
the molecular gas, we may find slightly lower maximum 
temperatures. However, this necessitates either the in- 
clusion of CO data at higher-J transitions or using so 
call ed molecular th ermometers (such as H2CO or NH3; 
see lAo et al.l l20l"ll ). High kinetic temperatures would 
be expected in several scenarios, among them shocks as 
well as heating through X-ray radiation from the AGN 
(e.g., iMeijering et al.l boOTi : iMeiiering fc SpaansI l2005l ). 
Nevertheless, a significant fraction (>>10-20%) of the 
warm molecular gas in this bridge regions seems to ex- 
hibit high kinetic temperatures which is a much larger 
fraction than so far found (for the overall molecular 
gas reservoir) in starbu rst and other Seyfert galaxies 
(<30%, see for instance iRoussel et al.l 120071 : iDale et all 
I20M iRigQPOulou et al.ll200l) . 

The density is also well constrained and is found to 
be in the range o f lO'^'^-lO^'^ cm~ '^ , consistent with 
previ ous findings (jKrips et al.l l2008l : IMatsushita et al.l 
I1998D . Considering the range of assumed abundance 
ratios, the column density of ^^CO approximately 
spans a range between ~10^^ ° cm~^ and 10^^ " cm~^. 
Even though CO might not be the best tracer o f 
the isotopic ratio (e.g., iMartm et all |2010[) 

we assume it to be a first approximation to this ra- 
tio. The carbon ratio of 26 found for the abso- 
lute lowest x^ is similar to the value of 20 measured 
in th e Galactic Center region (e.g., IWilson fc RoodI 
Il994[ ) and l ower than that derived towards nearby star- 
bursts (e.g..lHenkel et al.]ll993al : lHenkel fc Mauersberged 
Il993bt IHenkel et al.lll994D . Such enrichment would 
point towards a highly nuclear processing of the ISM in 
the central region of NGC 1068. 

5.3.2. HCN & HCO+ emission 



In IKrips et all (|2008D we carried out an LVG analy- 
sis for NGC 1068 based on the HCN and HCO+ single- 
dish lines ratios with kinetic temperatures not exceed- 
ing 200 K. As our new simulations indicate kinetic tem- 
peratures lying significantly above 200 K, we repeated 
the simulations with RADEX allowing for a larger range 
in kinetic temperatures. The line fluxes for the HCN 
and HCO+(J=1-0) emission are taken fr om PdBI obser- 
vatio ns at '^l" angular resolution (Garc fa-Burillo et al.l 
l2008i) which will be analysed in more detail in a later pa- 
per by Usero et al. (in prep.). The results of the RADEX 
simulations for HCN and HCO+ are shown in Fig. [22] 

Considering the restrictions for the kinetic tempera- 
tures (>200 K), we find solutions (i.e., with low x^) with 
RADEX for which the HCN-to-HCO+ abundance ratio 
lies between Xy^Q+~10-500 which is higher than that 

found in starforming/starbursting regions (X^^qI^—I). 
However, good solutions are also found for lower kinetic 
temperatures. This is due to the fact that HCN and 
HCO+ are not as sensitive to changes in the kinetic tem- 
peratures as '^^CO and ^'^CO. They are better indicators 
of changes in the volume density, as can be nicely seen 
in Fig.|22l the volume density is restricted by the x^-test 
to a very small area and independently yields values of 
n(H2)~103-5-104-5 cm-3 similar to the ^^CO and "CO 



10 



Krips et al. 



results. 

The simulations^^ indicate column densities for HCN 
in the range of N(HCN)~1012.o_;l013-5 cm'^ which is 
smaller (a factor of ~10) than what was found by 
iKrips etall (|2008l ) with the LVG code. However, results 
are quite similar to iKrips et al.l ()2008[ ) when assuming 
similar kinetic temperatures. 

Comparing the HCN column densities to those of ^^CO 
(N(i2CO)~10i7"-10i9 cm-2), we obtain abundance ra- 
tios between HCN and ^^CO of Xh^o>io3-5 ^hich seems 
to be still compatible with a slightly increased abundance 
of HCN. Comparing the column densities of HCO+ and 
^^CO, we find a somewhat decreased HCO"*" abundance 
(by a factor of at least 10 lower than found in galac- 
tic starforming r egions). This agrees well with previ- 
ous r esults (e.g.. iGarcia-Burillo et al.l l2010t IKrips et"ari 
[20081 lUsero et al.ll2004D suggesting an increased forma- 
tion (and hence increased abundance) of HCN due to an 
XDR in the center of NGC 1068. 

6. SUMMARY & CONCLUSIONS 

The SMA and PdBI observations of the (sub-)mm 
emission in NGC 1068 presented in this paper show a 
complex distribution, kinematics and excitation condi- 
tions of the molecular gas. The (sub)mm continuum and 
molecular line emission is interpreted as follows: 

1. ) The cm/mm-continuum emission seems to be best 

reproduced by electron-scattered synchrotron emis- 
sion. Thermal fre e-free emission as proposed by 
iHonig et al.l ()2008[) overpredicts the high angular 
resolution 1 mm continuum emission. 

2. ) The molecular gas is found to display a very com- 

plex kinematic behaviour in the ^^CO, HCN and 
HCO+ lines which is not reproducable by a tilted- 
ring model approximating a warped disk with cir- 
cular motions. Instead, a dominant rotating disk 
plus a radial outflow of some of the gas in the CND 

Pl ease note that we consider the two highest contours in 
Fig. 1221 as being acceptable solutions 



is proposed as an alternative explanation to ac- 
count for the non-circular motions. 

3.) The different line ratios from the ^^CO, ^^CQ, 
HCN, and HCO"'" emission seems to be consis- 
tent with moderately dense and warm gas, both 
being further support for a gas scenario in which 
heated and compressed by a shock (at least in the 
North/North-Eastern part of the ring). The high- 
est line ratios are found close to the AGN and/or 
jet-CND 'contact'-point. In this picture, the in- 
creased kinetic temperatures seem to be one of the 
culprits for the unusually high HCN-to-CO(J=l-0) 
line ratios due to a hypo-excitation of the C0(J=1- 
0) line emission. 

5.) Consistent with previous papers, we find further 
indications of an increased HCN abundance in 
NGC 1068 (by a factor of ~4-10), and a decreased 
HCO+ abundance (by a factor of '-^5-10), explain- 
ing the high HCN-to-HCO+ abundance ratio in the 
CND of this source. 
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TABLE 1 
Properties of NGC 1068 
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AGN type 

Dynamical Center . . 

Redshift 

Systemic Velocity . . 
Luminosity Distance 

Scale 

Inclination Angle . . . 
Position Angle 



NGC 1068 
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{R)SA(rs)b 
Seyfert 2 
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(1137±3)km s-i 
12.6 Mpc*" 
61 pc/arcsec 
40° 
278° 



NED" 

Khachikian fc Weedman ('1974 1 
Gallimore et al. (2001) 



Huchra et al. (1999) 
Huchra et al. (1999) 
NED" 
NED"''' 
Bland-Hawthorn et al. ('1997') 
Bland-Hawthorn et al. fl997 ) 



NED: NASA/IPAC Extragalactic Database 
^ Used cosmology: Hubble constant Ho = 73 km s~^/Mpc, Omega(matter) 
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TABLE 2 

Chronological Summary of observations carried out for NGC 1068 



Molecular 


Telescope 


Observing 


Config." 


Frequency 


Band 


Zenith 


Tsys 


RMS 


Synthesized Beam 


Line 




Dates 




at rest 




Opacity 




Noise'' 


majorxminor,P.A.'^ 






(YYYY-MM) 




(GHz) 




at 225 GHz 


(K) 


(mJy) 


("x",°) 


i3cO(J=2-l) 


SMA 


2008-01,2008-02 


EX 


220.399 


LSB 


0.10-0.20 


100-150 


12 


1.0x0.8,30'' 


& Ci80(J=2-l) 








219.560 


LSB 










HCO+{J=4-3) 


SMA 


2007-08,2007-09 


EX 


356.734 


USB 


0.06-0.13 


200-400 


32= 


1.0x0.8,30-^ 


i2cO(J=3-2) 


SMA 


2007-09 


EX 


345.796 


LSB 


0.06-0.13 


200-400 


51 


1.0x0.8,30-^ 


HCO+{J=3-2) 


SMA 


2006-11 


EX 


267.558 


LSB 


0.06-0.15 


100-150 


33 


1.0x0.8,30-'' 


HCN{J=3-2) 


SMA 


2006-01 


EX 


265.886 


LSB 


0.05-0.10 


100-150 


24<= 


1.0x0.8,30'='-'' 




SMA 


2006-10,2006-11 


VEX 




LSB 


0.10-0.22 


100-200 




0.53x0.46,30'*-'= 


"CO(J=3-2) 


SMA 


2005-10 


C 


356.734 


USB 


0.05-0.06 


200-300 


39 


2.4x2.1,28'* 


13CO(J=1-0) 


PdBI 


2003-02 


A 


110.201 


SSB 


n.a. 


150-300 


1.7 


2.5x1.9,28-'' 


i2CO{J=2-l) 


PdBI 


2003-02 


A 


230.538 


SSB 


n.a. 


200-600 


4 


1.0x0.8,30"= 



^ SMA configurations: C^compact (baselines up to 70 m), EX=extended (baselines up to 220 m), VEX=very extended (baselines 

up to 500 m); PdBI configurations: Almost extended (baselines up to 500m) 

^ in 17 km s~^ wide channels 

'= P. A. is measured from North to East 

using uniform weighting 
'= combined for all tracks 
' using natural weighting 



TABLE 3 

Continuum Parameters for NGC 1068. 



A 


Synth. Beam 


Ao" 


AS" 


Peak 


Flux 


Deconv. Size 




majorxminor,P.A. 






Flux 


Density'' 


majorxminor,P.A. 




{"x"°) 


(") 


(") 


(mjy/beam) 


(mJy) 


("x",°) 


1.4 mm 


1.0x0.8.30 


-l-0.13±0.02 


-l-0.25±0.02 


19±2 


28±3 


(0.9±0.1)x(0.7±0.1),(50±20) 


1.3 mm 


1.0x0.8,30 


-f0.18±0.02 


-l-0.23±0.02 


15±1 


22±2 


(0.6±0.1) X (0.5±0.1),(40±20) 


1.0 mm (NA)'= 


1.0x0.8,30 


-f0.17±0.05 


-l-0.15±0.10 


19±2 


24±3 


(0.8±0.1) X (0.4±0.1),(20±10) 


1.0 mm (UN)'* 


0.5x0.4,30 


+0.13±0.03 


-l-0.07±0.03 


12±2 


13±2 


(0.3±0.1)'= 


810 fim 


2.1x2.0,80 


-0.04±0.2 


-0.07±0.2 


30±4 


41±11 


(l.litO.l) X (0.8±0.2),(90±20) 


850 fj.m (NA) 


1.0x0.8,30 


+0.33±0.07 


-l-0.05±0.06 


24±2 


50±7 


(l.l±0.1)x(0.8±0.2),(90±20) 


850 fim (UN) 


0.6x0.5,30 


-f0.33±0.09 


-0.03±0.07 


16±3 


30±5 


(l.l±0.1)x(0.8±0.2),(90±20) 



" The offsets are with respect to aj20oo =02h42m40.70s and (5j20oo =-00°00'47'.'9 whic h is almost identical to the radio p osition 
of the AGN in NGC 1068 of aj2ooo =02h42m40.709s and 5j2ooo =-00°00'47'.'95 (e.g. iGallimore et al.l[200l : iKrips et al.l[2006 ). 
Positional errors are of pure statistical nature and were derived from the Gaussian fit to the data. They do not include absolute 
positional uncertainties from the calibration, which are estimated to be ~0'.'l. 

'' Flux errors are purely statistical and do not account for uncertainties of the flux calibration. The latter are estimated to be 
of the order 10-20% (see text). 

'= Averaged continuum emission derived from the HCN(J=3-2) (vex-(-ext) and HCO^(J=3-2) observations (ext). Data were 
mapped using natural weighting (NA). 

"* Averaged continuum emission derived from the HCN(J=3-2) observations alone (vex-|-ext) using uniform weighting (UN). 
*= Here, only a circular Gaussian fit has been carried out, while for the rest an elliptical Gaussian has been fitted to the data 
(see text). 
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TABLE 4 

Individual components of the molecular line emission in NGC 1068. 



Molecular 


Component 




AS"- 


Vel. Integrated 


Spatially Integrated 


Line 








Peak Intensity" 


Intensity" 






(") 


(") 


(Jy beam" km s""*^) 


fJv km s""*^") 


HCN(J=3-2) 


E-knot 


-|-1.0±0.1 


+0.1±0.1 


51±8 


110±20 


W-knot 


-0.8±0.2 


-f0.2±0.2 


32±4 


70±10 


HCO+{J=3-2) 


E-knot 


+0.9±0.1 


+0.0±0.1 


28±5 


52±5 


W-knot 


-0.8±0.2 


+0.5±0.2 


14±2 


40±5 


HCO+(J=4-3) 


E-knot 


+0.9±0.1 


+0.1±0.1 


27±6 


98±20 


i2C0(J=3-2) 


E-knot 


+1.1±0.1 


-O.litO.l 


470±20 


1330±200 


W-knot 


-1.2±0.1 


-O.litO.l 


270±10 


720±100 


i3CO(J=3-2) 


E-knot 


+1.1±0.1 


+0.3±0.1 


24±3 


40±10 


i2cO(J=2-l) 


E-knot 


+1.0±0.1 


-0.2±0.1 


70±10 


290±20 


W-knot 


-1.3±0.1 


-0.3±0.1 


30±5 


180±20 


"CO(J=2-l) 


E-knot 


+0.6±0.2 


+0.4±0.2 


10±2 


18±3 


W-knot 


-1.2±0.3 


-fO.litO.l 


8±2 


16±3 


Ci**0(J=2-l) 


E-knot 


+0.3±0.3 


+0.5±0.3 


7±1 


15±3 


13CO(J=1-0) 


E-knot 


-fl.0±0.4 


-0.4±0.4 


0.4±0.2 


0.5±0.2 


W-knot 


-1.3±0.4 


-0.8±0.5 


0.4±0.2 


0.5±0.2 


12CO(J=1-0) 


E-knot 


+0.5±0.1 


+0.1±0.1 


40ibl 


90±2 


W-knot 


-l.litO.l 


-0.5±0.1 


19±1 


40±2 



" The parameters were determined by fitting a one- or two-component elliptical Gaussian profile to the ut;-data of each line. 
Errors include the statistical uncertainties from the Gaussian fit and those from the calibration (~10-20%). Offsets are with 
respect to the center position specified in Table [S] 



TABLE 5 

Molecular line parameters derived from the different line spectra of NGC 1068. 



Molecular 


Velocity 


Line 


Line 


Vel. Integrated 


Vel. Integrated 


SD Beam 


Line 


OfTsct"'''''^ 


Flux" 


Width"'*'''* 


Intensity"''' 


SD Intensity'^ 


{") 




(km s-i) 


(Jy) 


(km s-i) 


(Jy km s-i) 


(Jy km s"-*-) 




HCN(J=3-2) 


-30±10 


0.63±0.08 


200±30 


150±10 


190±10 


9.5" 


HCO+{J=3-2) 


-40±10 


0.23±0.03 


190±50 


50±6 


80±8 


9.5" 


HCO+{J=4-3) 


-40±10 


0.24±0.03 


240±40 


60±6 


70±10 


14" 


i2C0(J=3-2) 


-30±3 


12.3±0.40 


170±10 


2130±5 


2600±300 


14" 


i3co(J=3-2) 


-40±10 


0.43±0.03 


230±30 


lOOilO 


170±20 


14" 


i2CO(J=2-l) 


-30±6 


2.20±0.10 


230±20 


529±2 


950±6 


12" 


i3CO(J=2-l) 


-lOitlO 


0.50±0.04 


60±10 


30±1 


55±7 


12" 


Ci«0(J=2-l) 


-l-3±10 


0.10±0.02 


50±10 


5.1±0.3 






12CO(J=1-0) 


-l-3±2 


0.50±0.02 


240±10 


120±4 


650±80 


21" 


i3cO(J=l-0) 


-l-8±10 


0.009±0.002 


140±30 


1.3±0.3 


56±7 


21" 



" The line parameters have been determined by fitting a single Gaussian line to the (spatially integrated) spectrum for each 
molecule. The line emission has been thereby integrated over the central 4" in NGC 1068. 

statistical error from the Gaussian fit only. 

with respect to vlsr=1137 km/s. 

Full Width at Half Maximum (FWHM) 

single dish (SD) integrated intens ities as measured with the I RAM 30m and t he JCMT telescope in the central 10-30" of 
NGC 1068 (taken from: |lsrael 2009i: iPerez-Beaupuits eralll2009l : iKrips et al.|[2008h . The values were converted to Jansky scale 
using S[Jy]/Tmb[K]=4.71 (30m) and S[Jy]/Tn,b[K]=15.6 (JCMT). 



TABLE 6 

Molecular line ratios for NGC 1068". (X[K]/Y[K]) 



Y 










13CO 


Ci80 


HCN 


HCO+ 




X 


J=l-0 


J=2-l 


J=3-2 


J=l-0 


J=2-l 


J=3-2 J=2-l 


J=l-0 J=2-l J=3-2 


J=l-0 J=3-2 


J=4-3 












E-knot 








i2co J=l-0 




0.3±0.2 


0.2±0.1 


20±10 












J=2-l 


2.9±0.3 




0.3±0.2 




5.0±3.0 


17±4 


4.0±2.0 


7.0±3.0 




J=3-2 


6.0±3.0 


4.0±2.0 








25±6 


5.0±3.0 


8.0±3.0 




J = i— U 


n n _Ln nu 
U.Uo±U.Uo 








0.2±0.1 


0.08±0.04 








T — 9 1 




n 9-1-0 1 
u. zmu. i 




O-U'i 
O.UlCo.U 




0.7±0.3 1.4±0.3 


'i-un 9 


f^-un 9 
■ ■ ■ u.omu.z 




T Q O 

J — o— z 






O 0/1 -UO 01 


1 'i-\-A 


1.4±0.4 




o 9-un 1 


o c;-un 'i 

■ ■ ■ U.OltU.vS 




/^18/^ T O 1 

yj J — z— 1 




n ntj-un riQ 
U.UDitU. Uo 






0.8±0.3 










IlLjiN J — i-—\J 




















J — o— z 




U.oitU.l 


fl 1 /I -UO 0*3 




4±1 


4.0±1.0 








rlL/U ^ J — i— U 




















J — o— Z 




n 1 Q-|_n 09 


nv-uo ni 
u.u 1 mu.ui 




2.0±0.4 


2.0±0.5 


■ ■ ■ ■ ■ ■ u.omu.o 






















■ ■ ■ u.omu.o 














W-knot 








"CO J=l-0 




0.8±0.7 


0.3±0.2 


50±30 












J=2-l 


1.5ib0.5 




0.3±0.2 




10±2 




3.4±0.6 


e.Oiti.o 




J=3-2 


4.0ib2.0 


3.0±2.0 








25±6 


e.Oiti.o 


11±2 




13 r^/^ Tin 


n no 1 n ni 








0.13±0.06 










T O 1 

J — z— i 




n 1 0-Un 09 




fi O-U/l o 






o /1-un 9 

■ ■ - ■ ■ - U.4itU.Z 


■ ■ ■ U.DltU.o 




T Q O 

J — o— Z 






o O/i -UO m 














O U J=z— i 




















wnivr T 1 n 

rn^iN J — i— u 




















T Q O 

J — o— z 




n Q-Un 1 


n o_un 1 




3.0±2 






1 n_l_n K 




XJf~*(~\-\- Tin 




















T Q o 

J — o— z 




O 1 7-I-0 OQ 


n oo-Uo no 
U.uyitU.Uz 




1.7±0.4 




■ ■ ■ ■ ■ ■ U.OltU.O 






T A r> 

J — 4— o 






























Total 








i2co J=l-0 




0.3±0.2 


0.3±0.2 


40±20 












J=2-l 


2.0±1.0 




0.3±0.2 




8±3 


94±6 


4.0±2.0 


9.0±4.0 




J=3-2 


4.0±2 


4.0±2.0 








41±2.0 


8.0±4.0 


16±1 




i^CO J=l-0 


0.04±0.02 








0.17±0.04 


0.2±0.1 








J=2-l 




0.2±0.1 




e.Oiti.o 




1.0±0.5 


0.3±0.2 


0.6±0.3 




J=3-2 






0.025±0.001 


4.0±1.0 


0.7±0.3 




0.2±0.1 


0.4±0.2 




Ci^O J=2-l 




O.OlliO.OOl 






0.2±0.1 










HCN J=l-0 




















J=2-l 




















J=3-2 




0.3±0.1 


0.12±0.01 




3.0±0.3 


5.0±3.0 




1.9±0.2 




HCO+ J=l-0 




















J=3-2 




O.llitO.Ol 


0.08±0.02 




1.8±0.1 


3.0±0.2 


0.5±0.3 




1.6±0.8 


J=4-3 












1.6±0.1 




0.6±0.3 





The lino ratios were derived by spatially averaging over the respective region from the velocity intergated line ratio maps from Fig. 1121 The errors denote thereby the 
standard deviation from the averaged values. Please note, that in some cases, the line ratios might vary from those estimated from Tables[4]and[5] See text for a discussion. 
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Rodio position of AGN: white cross (G04, VLA) _ _ Rodio position of AGN: winite cross (G04, VLA) 

_ 1.3mm cent.: grey scale & contours (PdBI) 1.3mm cont.: grey scale & contours (PdBI) . 

1.4mm cont.: black contours (SMA 1.0mm cent.: black contours (SMA; NA) 

1.0mm cont.: white contours (SWA; UN) 



:(c) 





:(b) 



Rodio position of AGN: white cross (G04, VLA) J 
1.3mm cont.; grey scale Sc contours (PdBI; 
850fim cont.: black contours (SMA, NA) 
850/j.m cont.: white contours (SWA, UN) 





- 0.005 





RA - offset (") 





offset (") 





RA - offset (") 



Fig. 1. — Continuum emission of NGC 1068 at A=1.4 mm (Fig. la; black contours), 1.3 mm (Fig. la-c; grey scale and grey contours), 
1.0 mm (Fig. lb; black contours) and 850 fj,m (Fig. Ic; black contour s), observed with the SMA and the IRAM PdBI. The white cross 
denotes the position of the AGN measured by iGallimore et al.l II2004I . G04). The contours of the 1.3 mm continuum emission (PdBI) 
start at 5(t=4 mjy in steps of Icr. a) The contours of the 1.4 mm continuum emission (SMA) start at 5cr=4 mjy in steps of la. b) The 
contours of the 1.0 mm continuum emission (SMA, NA) start at 3(t=1.6 mjy in steps of Icr, while the contours of the uniformaly mapped 
1.0 mm continuum emission run from 3(t=2.3 mjy in steps of Itr. c) The contours of the 850 fira continuum emission (SMA, NA) start at 
3(7=2.4 mjy in steps of Icr, while the contours of the uniformaly mapped 850 /xm continuum emission run from 3cr=2.6 mjy in steps of Icr. 
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~i — ' — ' — ' — \ — ' — ' — ' — \ — ' — ' — ' — I — ' — ' — ' — r 
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AGN radio po&itinbri'sjr^oss'- (S^iA/¥tA) 



— I — I — I — I — I — I — I — I — I — I — I — I — I — 1— 
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jO(2-1); black contours (PdBI) 
GN radio position: cross (SMA/VLA) 



— I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1- 
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'^C0{3-2); black contours (SMA) 
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__ 

^ r- ■ ■ " ■: 



(SMA) 



a- 




'^C0(2-l): blflcrT^ntours (SMA) 

OSS (SMA/VLA) 

o 




400; 




P t ' ' •' — I — ^ 
-GQp— .2:-): gr^-"SEalB' + co^^^., ir-= rsw. 
C"0(2-i;; blact' CLir-itoura (SMA) 
A6fJ radio po^-ttiijn; ^roee i SMA/VLA} 



-t- 



-C0(3-?.) c^=^\:^7r5l^: i- con to I ITS {"5 MAT 
'^CO(I-O): black contours (PdBI)' 
AGN radio position: cross (5MA/VLA) 

C) 



[^W-O): black contours (PdBI) 
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n. 0- 
I 



-4 




X9> 



•(hi 




40o: 




S: 



f?P ~ uffset (■') 



:0' —Z, 
Prt - offiSt t") 



,9, 



Fig. 2.— Velocity integrated line emission of l2cO(J=3-2) (a-g), HCN(,I=3-2) (a), HCO+(J=3-2) (6), HCO+(J=4-3) (c), l2cO(J=2-l) 
(d), i-^CO(J=2-l) (e), i3CO(J=3-2) (/), Ci**0(J=2-l) (</), i2co(J=l-0) (ft) and i«CO(J=l-0) (j) in NGC 1068, observed with the SMA 
and the IRAM PdBI. Contour levies are; "co(J=3-2) - from IOct by 6cr with lo-=4.8 Jy km s^^; HCN(J=3-2) - from Sct by la with 
lcr=2.6 Jy km s^^; HCO+(J=3-2) - from 2cr by Icr with 1(t=3.4 Jy km s'^; HCO+(J=4-3) - from 3a by Icr with la=2A Jy km s^i; 
12C0(J=2-1) - from 5a by 5cr with lo-=1.2 Jy km s^^; i"'CO(.]=2-l) - from 3cr by Icr with lcr=0.9 Jy km s^^; i-^CO(J=3-2) - from 5a by 
1<T with 1(T=3.5 Jy km s"!; Ci*0(J=2-l) - from 3a by la with 1<t=0.9 Jy km s"!; i2cO(J=l-0) - from 5a by 3a with 1(t=0.4 Jy km s'^; 
"CO(J=1-0) - from la by la with 1ct=0.1 Jy km 8"^ 
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0.4 



0.2 




0.4 - 



0.2 



Av (km/s) 



-200 

Av (km/s) 



200 



Fig. 3. — Spatially integrated spectrum of different molecular lines in NGC 1068. The single-line Gaussian fit is indicated with a dotted 
red line (parameters are listed in Table [Sjl while the multiple Gaussian fit is plotted with a dashed blue line. 
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C0(2-1) 

tota 
E-knot 
W-knot 



^ I ' 
NGC1068 



HCN(3-2) 

— tota 

■■■ E-knot 

- W-knot 



HC0''(4-3) 

- tota 

■•■ E-knot 

- W-knot 



IJ 



NGC1068 




(e) 



NGC1068 



(f) 



-200 200 -200 200 

Av (knn/s) Av (km/s) 

Fig. 4.— Spatially integrated spectrum of the 12CO(J=1-0), (talccn from ISchinnerer eraL|[2000h . i2CO(J=2-l) and ^^CO{J=3-2) {left 
column) and HCN(J=3-2) and HCO^(J=4-3) emission over the E-knot {dotted blue) and W-knot {solid red) component of NGC 1068. 
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I 1 , 1 I 1 , 1 I , 1 , I , 1 r 

.v^so— velocities: red + blue cbntours (5MA) 
_^,mm continuum: cross (SMA/PclB!) 



i 




(a) 



1^— velocities; red + blue contours (5MA) 
rrim continuum: cross (SMA/PdBI) 



■(b) 

_J L 




1 , , ^ ,— 1 , I , , , I , , r 

- C0!7-!i: artH- scoift + conto^"-' ;"MA) 
so— velocities: red + blue contours (SMA) 
1mm continuum: cro^s (SMA/PdBI) 




Iso-velocWes; red + blue contours (StviA) 
'mm continuum: cross (SMA/PdBI) 



Iso— velocities: red + blue contours (SMA) 
1 mm continuum: cross (SMA/PdBI) 





Iso— velocities: red + blue contours (SMA) 
1 mm continuum; cross (SMA/PdBI) 





(d) 
4 



2 0-2 
RA - offset (") 
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Fig. 5. — Iso-velocity maps of the different molecular lines observed in NGC 1068. The grey scale correspond to the velocity integrated 
emission of each line with the same contours as used in Fig. [2] The velocity contours are in steps of 10 k m s~ ^ around the systemic velocity 
of NGC 1068. The grey lines indicate the cuts along which the position-velocity diagrams (see Fig. IIOI I were taken for the respective 
molecules (12C0(J=2-1),12C0(J=2-1) & HCN(J=3-2)). 
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Fig. 6. — Red- and blueshifted emission of NGC1068. Contours are in steps of a) from 7a by 5cr with lcr=2.0 Jy km s~^; b) from 5(T by 
5(T with lo'=0.9 Jy km s~^; c) from 4a by la with 1(7=0.5 Jy km s~^; d) from 3a by la with 1(t=1.9 Jy km s~^; e) from 3a by la with 
1(7=2.3 Jy km s~^; f) from 2a by la with 1(7=1.4 Jy km s~^. 
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-00°00'45" - 




Fig. 7. — Channel maps of the i2co{J=2-l) emission in NGC 1068. Please note that the i2cO(J=2-l) data were resampled to match 
the spectral resolution of the ^■^CO(J=3-2) data and facilitate a comparison, especially with respect to the line-ratio channel map in 
Fig. 1131 Contour spacing is in steps of 5(t=37.4 mjy beam~^. We use a spectral resolution of ~7 km/s. The zero channel corresponds to 
VLSR=1137 km of NGC 1068. 
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Fig. 8. - Channel maps of the i2cO(J=3-2) emission in NGC 1068. We use a spectral resolution of ~7 km s ^ and the original spatial 
resolution from the observations. Contour spacing is in steps of 5<t=400 mjy beam"'^. The zero channel corresponds to vlsr=1137 km s~^ 
of NGC 1068. 
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Fig. 9. — Channel maps of the ^''CO{J=2-l) emission in NGC 1068. Wc use spectral resolution of ~7 km s ^. Contour spacing is in 
steps of 2(T=11 mjy beam~^. The zero channel corresponds to VLgj^=1137 km s^'^ of NGC 1068. 
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Fig. 10.— Position-velocity diagram of NGC 1068 for the i2CO(J=2-l) (grey scale), i2CO(J=3-2) (black contours, top figure) and 
HCN(J=3-2) (black contours; bottom figure) emission. 
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Fig. 11.— Parametrisation of the kinematics in the CND of NGC 1068 for the i2cO(J=2-l) emission (filled black data points) taken 
along the slits at different position angles (at 0, 30, 60, 90, 120, and 150°). cr is the average line-of-sight velocity dispersion, e the rms 
variation of the velocity from point to point and A the rotational velocity along each slit (see text in Sectior i5T2l for a more detailed 
description). VNRot = violent non-rotators, CNRot = calm non-rotators and Rot = rotators. The filled grey data points were determined 
from the model discussed in Section 15.21 
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Fig. 12. — Velocity integrated molecular line intensity ratios for NGC 1068 above a >lcr threshold. The white cross marks the po- 
sition of the mm continuum emission that is associated with the AGN. The different line ratios are: a) ^^CO(J=3-2)-to-^-^CO(J=l-0); 
b) i2cO(J=3-2)-to-i2cO(J=2-l), c) i2cO{J=2-l)-to-i2cO(J=l-0), d) i2cO(J=3-2)-to-i3cO{J=3-2), e) HCN(J=3-2)-to-i2cO(J=3- 
2), /; HCO+(J=3-2)-to-i2cO(J=3-2), g) HCN(J=3-2)-to-HCO+{J=3-2), h) HCO+{J=3-2)-to-HCO+ {J=4-3), i) i3CO(J=2-l)-to- 
i2CO(J=2-l), j) i3cO{J=2-l)-to-HCN(J=3-2), k) i3cO{J=2-l)-to-HCO+ (J=3-2), I) i3cO{J=2-l)-to-Ci80{J=2-l), m) i3cO(J=l- 
0)-to-i3cO(J=2-l), n) i3cO(J=l-0)-to-"cO(J=l-0), o) "cO(J=2-l)-to-"cO(J=3-2), p) "cO(J=l-0)-to-i2cO(J=l-0) 
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Fig. 13.— i2CO(J=3-2)-to-i2co(J=2-l) line ratio above a 4(7-threshold for each line. 
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Fig. 14. — Spatially averaged line ratio for the eastern and western knot as function of velocity, derived from Fig. llijI The error bars denote 
the variance of each averaged value, which does not exceed 50% in most cases. The missing values for velocities —120 and —200 km s~^ 
and +100 to +140 km s~^ are due to the lack of emisison in the respective knot. The dashed lines represent the median of the line ratios 
for the two knots. 
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+ (Sub-)mm data of NGC1068 (this paper) 
: ®IR/cm data of NGC'068 {Haenig et al. j008; 
- O Data observed with angular resolution of <r' 




CSub-)mm data of NGC1068 (this paper) 
®IR/cm data of NGC1068 (Hoenig et al, 2008) 
O Data observed with angular resolution of <!' 




(Sub— )mm dato of M -i -hio paper) 

■V IR/cm data of Nf- - .1 . 

O Data observed with orgular resolution of <1' 




Fig. 15. — Spectral ener g y dist ribution of the continuum emission in NGC 1068, based on data from this paper. iKrips et ah! l|2006l ') (black 
crosses), and lHonig et all l|2008l ) (grey crosses). The dotted line respresents the model for the radio continuum emission (either electron 
scattered synchrotron emission (a), synchrotron emission (&), and free-free absorption (c)), the dashed line represents the model for the 
IR data (a two-temperature grey body a-c), and the solid line represents the composite of both (a-c). The data observed at an angular 
resolution below 1" are additionally marked with a circle. 
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Fig. 16. — Velocity channel maps of the CO model compared to the l2cO(2-l) emission. 



Millimeter observations of NGC 1068 



31 




32 



Krips et al. 



2 - 



CO 

^ 



o 
in 

Q 



- Observations 



4 - 



2 - 



o 

Q 



4 - 



LlI 



- 



Moment — 



Moment- 



- Observations 

H — I — I — I — I — ^ 
Moment — 2 



n ' ' ' \ ' ' ' \ ' \ ' \ ' ' ' r 

C0(2-1) Moment-0 



co(; 





- fvlodel 



H — I — h 



H — I — h 



C0(2-1) Moment-1 




- Model 
^ — \ — \ — \ — \ \ — \ — I — I — I — I 

00(2-1) Moment-2 




Observations 



Model 



_L 



_L 



2 0-2 
RA - offset (") 



-4 



H — h 



00(2-1 ; 




I — I — I — \ — 

00(2-0 




i -2 

RA - offset (") 



-4 



Fig. 18. — Moment maps of the CO model compared to the ^^CO(2-l) emission. The velocities are plotted in steps of 20km s ^ for the 
Moment-1 and Moment-2 maps. 
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Fig. 19. — Blue- and Redshifted emission of the CO model compared to the ^^CO(2-l) emission. 
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Position-velocity diagram of the CO model compared to the ^^00(2-1) emission. 
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Fig. 21. — x^-&t results obtained from the RADEX simulations of the excitation conditions of the molecular gas. Shown are four different 
[^■^CO]/[^^CO] abundance ratios (=10,26,52,110) for three different ^^CO column densities respectively. The middle panel shows the best 
X^-fit found for each abundance ratio. 
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Fig. 22. - x^"fit results obtained from the RADEX simulations of the excitation conditions of the molecular gas. Shown are four different 
[IfCN]/[HCO+] abundance ratios around the standard galactic value of [IfCN]/[HCO+]~10 for three different HON column densities 
respectively. The middle panel shows the best x^-^^ found for each abundance ratio. 



